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ABSTRACT 

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are 
now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us 
from the surface of last scattering, the CMB photons are expected to undergo modifications induced by 
the intervening large-scale structure. Among the expected secondary effects is the weak gravitational 
lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that 
uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the 
power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions 
from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the 
first direct constraints of this lensing potential power spectrum and find that it has an amplitude of 
A L = 0.96 ± 0.60, 1.06 ± 0.69 and 0.97 ± 0.47 using the W, V and W+V bands, respectively. 
Subject headings: cosmology: cosmic microwave background — cosmology: observations — cosmology: 
theory — gravitational lensing 



1. INTRODUCTION 

Measurements of Cosmic Microwave Background 
(CMB) anisotropics have served as the strongest exper- 
imental probe of the early Universe to date (Spergel et 
al. 2003). Temperature fluctuations capture the physics 
of the primordial photon-baryon fluid undergoing oscil- 
lations in potential wells sourced by primordial density 
perturbations (Hu et al. 1997; Hu & Dodelson 2002). 
The anisotropics in the CMB also contain information 
related to the Universe at late times, as CMB photons in 
transit to us encounter large-scale structure. The main 
mechanisms that affect the frequency and direction of 
propagation of CMB photons are the gravitational inter- 
action with time varying potential perturbations (Sachs 
& Wolfe 1967; Rees & Sciama 1968), and Compton scat- 
tering by electrons due to the reionization of the uni- 
verse (Sunyaev & Zel'dovich 1970). These mechanisms 
create specific signatures in the CMB that can be used 
to extract properties of the large-scale structure. 

One observational signature present in the CMB is the 
gravitational lensing modification by the projected dark 
matter distribution integrated along the line of sight (Hu 
2000; Lewis & Challinor 2006). Thi is analogous to 
weak lensing of galaxy shapes, now understood to be a 
strong probe of the dark matter distribution (Wittman 
et al. 2000; Bacon et al. 2000, van Waerbeke et al. 
2000; Bartelmann & Schneider 2001). Unlike the case 
of lensing measurements with galaxy shapes, which re- 
strict studies out to the era of z ~ 2, lensing of the CMB 
traces all the way back to the surface of last scatter- 
ing at z ~ 1100. The dominant contributions to CMB 
lensing arise from z ~ 2, but with a 30% contribution 
at higher redshifts (Lewis & Challinor 2006). From the 
ability to probe distance ratios and the integrated mat- 



ter power spectrum at early times, CMB lensing is un- 
derstood to be a powerful probe of certain cosmological 
parameters such as the neutrino mass and early dark 
energy (Kaplinghat et al. 2003, Smith et al. 2008, Les- 
gourgues et al. 2006, Li & Cooray 2006). 

Unlike lensing measurements with discrete galaxies, 
the primordial CMB sky is a continuous field, and a dif- 
ferent technique is therefore needed. Due to lensing de- 
flections, the temperature anisotropy O(n) = AT/T(n) 
measured on the sky becomes 9(n+a), where a = V</>(n) 
is the deflection angle due to lensing given by the angular 
gradient of the lensing potential. The lensing potential 
is the line of sight projection between us and the last 
scattering surface of the gravitational potential nr) 
in the Universe: 

</>(n) = -2 r dr dA{r ° ~ r) M r, nr), (1) 
Jo d A (r)d A (r ) 

where d A is the comoving angular diameter distance and 
r is the comoving conformal distance from the observer. 
Using this expression, we can consider the dominant cor- 
rection to the CMB temperature as a term that couples 
the deflection angle V</>(n) to the gradient of the tem- 
perature: 

6(n) = 0(n + a) ~ O(n) + V 4 0(n)V l O(n), (2) 

where 0(n) is the unlcnscd temperature map. 

The CMB lensing effect is a modification to the an- 
gular gradient of temperature on the sky. If the CMB 
were completely isotropic, lensing modifications would 
not leave a change as lensing conserves surface bright- 
ness. Further, the lensing signatures are at second-order 
in temperature leading to a distinct non-Gaussianity pat- 
tern at the four-point function level, or trispectrum in 
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Fourier space (Zaldarriaga 2000; Hu 2001b). Quadratic 
statistics can be devised to probe the gradient structure 
of the CMB temperature and to extract the projected 
line of sight dark matter density field (Hu 2001a; Cooray 
& Kesden 2003; Okamoto & Hu 2003). Attempts have 
been made for an indirect signature of the lensing ef- 
fect in all-sky maps by Wilkinson Microwave Anisotropy 
Probe (WMAP) three-year data by estimating the de- 
flection field and then cross-correlating that with a fore- 
ground density tracer field, such as NRAO VLA Sky 
Survey (NVSS) and Sloan Digital Sky Survey (SDSS). 
The overall signal-to-noise ratio for the measurement is 
2.5a (Hirata et al. 2008) to 3.4cr (Smith et al. 2007). 

With observations over seven years now completed, 
WMAP maps have improved in sensitivity to the extent 
that a direct measurement of the lensing signal can be 
pursued. We first construct an optimized statistic that 
probes the non-Gaussianity pattern induced by lensing 
via the trispectrum of the CMB temperature. We mea- 
sure this by four maps weighted differently and taking a 
power spectrum of the temperature squared field. In ad- 
dition to the lensing signal, this power spectrum contains 
a term associated with the Gaussian sky. We combine 
the constraints in data with a large suite of simulations, 
involving both the Gaussian temperature maps without 
lensing as well as maps with lensing included. 

This Letter is organized as follows. In the next Section, 
we outline the new estimator and how to weight CMB 
maps to extract the lensing signal. In Section [3] we show 
that our estimator works correctly for simulated data. In 

Section^we describe our analysis used to constrain Cf^, 
the detection of the lensing amplitude and describe null 
tests used as a sanity checks. In Section[5]we present our 
results and conclude with a discussion on implications 
for the future. 

2. DERIVATION OF ESTIMATOR 

We derive our estimator using the trispectrum of the 
CMB. (Hu 2001b; Okamoto & Hu 2003; Smidt et al. 
2010). In harmonic space, the lensing perturbed temper- 
ature field becomes (Eq. [2j: 
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x~[L(L + l) + l'(l' + l)-l(l + l)], 

where the term in large parenthesis is the Wigner-3j sym- 
bol. 

From here, it is straightforward to show that lensing 
does not generate a three point correlation function or 
a bispectrum. The first correction to non-Gaussianity 
from CMB lensing appears in the four point correlation 



function or, in Fourier space, the angular trispectrum: 
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where the first term on the right hand side represents 
the contribution from a Gaussian sky and is subtracted 
out in our analysis (Hu 2001b, Smidt et al. 2010). The 
second term is called the connected piece containing the 
contribution to the trispectrum from lensing and is ex- 
panded as 
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Here P/ 1 / 2 (L) is known as a "pairing" of the full trispec- 
trum (Hu 2001b). 



From this, we can write the full trispectrum ^ (L) 
keeping all the permutations as: 
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where the curly bracket is the Wigner 6-j symbol. 

It has been shown that for weak lensing of the CMB, 
the full trispectrum is approximated well by its pair- 
ing and thus we will approximate the full trispectrum 



^ TlllHL) 



P^t(L) (Hu 2001b). This simplifies our 



analysis significantly and we have verified that this ap- 
proximation is accurate to better than 3%, which is ad- 
equate given our overall constraints are only accurate to 
25%. 

To extract Cf from CMB maps, we develop a 
quadratic estimator based on the kurtosis power spec- 
trum of the CMB. Starting with Eq. [7j we find the 
trispectrum is of the form: 
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This lets us break up T^(l) into 4 x 9 = 36 distinct 
pieces: 
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allowing us to write K, l ' as 
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Here C; = Citif + Ni where C; is the temperature power 
spectrum, bi is the beam transfer function related to the 
experiment, Ni is the noise power spectrum is the tem- 
perature power spectrum and 1 l 2 (Z) is the full trispec- 



trum. Also, (I) is the full trispectrum with Cf 9 = 1 
leaving the righthand side proportional to Cf^. For con- 
venience we introduce the notation /C ; = Cf^ 'Kf . 

Each term of Eq. [nj f[^ 3l4 (l), can be broken up into 
the form 

f^(L) = h llhL h hUL F L a^ti?^5^. (13) 
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This allows us to recover JC\ ' from weighted CMB 
maps analogously to what was done for primordial non- 
Gaussianity by Munshi et al. (2009) and Smidt et al. 
(2010). 

First we introduce the necessary weighted CMB maps: 
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Here, the specific aj , /3 ; , 7 ; and <5 ; weightings for 
each of t he 3 6 terms of the trispectrum may be deduced 
from Eq. [13] and ai m comes from the CMB map being 
weighted. 

From these maps we can form a quadratic statistic for 
each term 



j^(x)(2,2) _ j^A ( ^ BW,G<") _D< x » 



(18) 



(21 



^(x)^(x) 



q(x)j~)(x) 



Im 



where [iwflf 1 )]^ are the spherical harmonics from the 
product of the and fiW maps with [G^D^] im 

defined similarly. 
From these definitions it is straightforward to recover: 
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Fig. 1. — The dashed line is the analytical calculation kurtosis 
(2 2) 

power spectrum ' ' for the fiducial cosmology in the WMAP 
W-band. The measured values are extracted from simulated lensed 
maps generated by Lenspix seeded by Cf . 

so that using Eq. [TT] we indeed recover Eq. [12] noting 
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In this form, the information from CMB maps can been 

extracted and Cf^ can be fit using Cf^ = k\ 2 ' 2 ^ //c| 2,2 - > . 
The signal-to-noise ratio for this estimator may be cal- 
culated as v/Ei(2Z + l)/C, (2,2) /24 where the factor of 24 
is needed to avoid over-counting of modes. We find the 
total signal-to-noise for each WMAP channel is 1.61 and 
when the two channels are combined is 2.5. 

3. TESTING THE ESTIMATOR 

To test this estimator we create 400 lensed and un- 
lensed CMB temperature maps for analysis. To create 
the unlensed maps, we calculate a temperature power 
spectrum from CAMB (Lewis 2000) from the best-fit 
WMAP 7- year parameters (Komatsu 2010). We then 
use the synfast Healpix (Gorski et al. 2005) routine to 
create a Gaussian map from this power spectrum and 
repeat until 400 maps are generated. These maps retain 
information from the power spectrum up to I — 1000 and 
have n s idc = 512. The anafast Healpix routine masks 
these maps with the KQ75 mask and produces the aim's 
to I = 750. We will refer to these a; m 's now collectively 

In addition, we create 400 noise maps for each each of 
our frequency bands V and W. We generate these maps 
from white noise as described in (Smidt et al. 2009) 
and denote the spherical harmonic coefficients as af m . 
The full simulated maps combine the Gaussian maps and 
noise maps in harmonic space by af m = a^Lbi + af m , 
where bi was obtained from the WMAP Teaix H and af m 
represent the spherical harmonics for the fultunlensed 
simulation. 

Next, we generate 400 Lensed CMB maps using Len- 
spix (Lewis 2008) that first creates Gaussian maps from 
an input temperature spectrum, then lenses the map 

1 http : //lambda. gsfc.nasa.gov/product/map/current 
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with an input lensing power spectrum. For both spec- 
tra, we used CAMB generated fiducial temperature and 
lensing power spectra based on the WMAP 7-year pa- 
rameters. From our maps, we obtain the lensed spherical 
harmonics a\ m using anafast. We then compile the full 
lensed simulation by adding these with the beam trans- 
fer functions and noise multipole moments as described 
above, utilizing the KQ75 mask. 
After the Gaussian and lensed simulations are ob- 

(2 2) 

tained, K.) ' may be extracted. Each simulation is 
weighted in the several combinations from which the kur- 



tosis power spectra IC 



(x)(2,2) 



are formed and then added 
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together to obtain the full estimator as shown in Eq 
These power spectra initially contain the Gaussian piece 
that must be subtracted out. Therefore, we take the 
mean of the kurtosis spectra obtained from the Gaussian 
simulations and subtract this from each kurtosis spec- 
trum obtained from the simulated lensed maps. To cor- 
rect for the cut sky, we use the technique developed by 
Hivon et al. (2001) that removes masking effects using 
a mode-coupling matrix akin to what was done in Smidt 
et al. (2010). The leftover connected pieces are averaged 

over yielding the theoretical K,f^ for fiducial cosmology. 

The kurtosis power spectrum Kt may be computed 



analytically from Eq. 9]|T2 using an unlensed temperature 



power spectrum and (7, obtained from CAMB. A plot 
of IC[ 2 ' 2 ^ versus C\ 



K^ 2 ' 2 ^ is given in Figure 1 der 
strating the estimator works well with simulated data 



To extract IC 



%2) 



4. ANALYSIS 
from data we use the raw WMAP 7- 



year Stokes-I sky maps for the V (60.8 GHz) and W (93.5 
GHz) frequency availabl^jl from the LAMBDA website. 
We analyze these maps with anafast (KQ75 mask) to 
generate the multipole moments for each frequency band 
out to l max — 750. We make no separate correction for 
unresolved point sources as their contamination for lens- 
ing measurements has been shown to be negligible (Smith 
et al. 2008). The correction for the cut sky is handled 
as described above. , s 

(2 2) 

These data maps are then weighted and K\ ' ' 
tracted with the Gaussian piece removed.. With K.) 

Cf+ = K 



(2,2) , ^(2,2) 



IK 



is ex- 

-(2,2) 

) 

can now be constrained. To con- 
strain each bin, we minimize \ 2 making use of the co- 
variance between bands and I bins as described in Smidt 
et al. (2009). The covariance matrix was computed and 
used to capture correlations between each Z-mode. A plot 
of these constraints is given in Fig|| For this fit, K, (2 - 2) 
is binned in I with 51 = 150. 
To avoid a false detection, we perform a null test 

measuring JC^ 2 ' 2 ^ on the difference of temperature maps. 
Here, we extract our estimator from the difference of the 
V and W frequency band maps. The reason for calcu- 
lating this difference is that the lensing signal should be 
the same in both maps, and therefore a subtraction of 
the maps should remove any lensing signal present. The 
difference is performed both on the raw Stokes-I WMAP 
7-year data maps, as well as the Gaussian simulations de- 
scribed above. Figure 2]inset shows no lensing is detected 
from these differenceclrnaps. We also tested our estima- 
tor against unlensed maps seeded with primordial non- 




FlG. 2. — The constraints on Cf using the K,( 2 ' 2 ^ estimator. In 
the larger plot, the red squares are the constraints by the V-band 
and the green circles are for the W-Band. The result of the null 
test is given in the smaller plot. These error bars represent la 
errors. 



Params. 


WMAP7 


WMAP7+A L 


WMAP7+ A L +Cf' p 


w 3 n b h 2 


22.51 ± 0.62 


22.59 ±0.63 


22.60 ±0.58 


10 2 fi DM h 2 


11.08 ±0.57 


11.04 ±0.54 


11.09 ±0.54 


T 


0.089 ± 0.016 


0.090 ± 0.015 


0.089 ±0.015 


n s 


0.967 ±0.015 


0.968 ±0.014 


0.968 ± 0.014 


n A 


0.734 ±0.031 


0.737 ± 0.028 


0.735 ±0.027 


Age/Gyr 


13.8 ±0.14 


13.7 ± 0.14 


13.7 ± 0.13 


Hi 


71.0 ± 2.7 


71.3 ±2.5 


71.1 ±2.4 


A L 


1.0 


0.87 ± 1.05 


0.97 ± 0.47 



TABLE 1 

Constraints on the cosmological parameters using 
CosmoMC. From left moving right: WMAP 7-year data 

ONLY; WMAP7 WITH A L ALLOWED TO RUN; WMAP7 WITH THE 
CONSTRAINTS ON Cf^ COMING FROM OUR /Cj 2 ' 2 ' ESTIMATOR FROM 
THE COMBINED W-BAND AND V-BAND DATA. THE 1 UNITS ON Ho 
ARE (KM S _1 MPC -1 ). 



Gaussinainity of / NL = 100 ( Eisner & Wandelt 2009) 
and found these are consistent with the Gaussian case 
for our estimator. This agrees with calculations from 
Lesgourgues et al. 2005. 

With constraints on Cf^, we use CosmoMC (Lewis 
and Bridle 2002) to constrain the lensing amplitude Al, 
as compared to that from the WMAP7 measurements of 
the temperature and E-mode polarization spectra C ; ee , 
Cf E , and Cf E alone. The parameter Al can be thought 
of as a measure of the degree of lensing in the CMB, 
where Al = 1 represents a universe with the expected 
amount of lensing signal and Al = is the case with 
no lensing. We sample the following seven-dimensional 
set of cosmological parameters: the baryon and cold dark 
matter densities Vt^h 2 and floAih 2 , the ratio of the sound 
horizon to the angular diameter distance at the decou- 
pling, 6 S , the scalar spectral index n s , the overall nor- 
malization of the spectrum A s at k = 0.002 Mpc -1 , and 
the optical depth to reionization, r. 

5. RESULTS AND DISCUSSION 
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Fig. [2] shows the scale dependent constraints on the 

(22). 

lensing power spectrum obtained from the K,\ ' estima- 
tor. We find the lensing signal to be compatible with 
the fiducial expectation. To our knowledge, this result 
provides the first direct constraints on Cf using CMB 
data encoding information for the matter distribution of 
the universe back to z ~ 1100. Furthermore, this mea- 
surement does not appear to be biased by instrumental 
effects since the null test described above, is compatible 
with zero. 

The results for measuring the lensing amplitude Al 
using the CosmoMC analysis described in Section [4] are 
summarized in Table [4j We first reproduce values for 
the cosmological parameters consistent with Komatsu et 
al. (2010) for the case with no constraints on Cj and 
Ai = 1. Next we carry out the same run, without a 
constraint on Cj , this time allowing Al to vary and 
find that WMAP 7-year data alone is consistent with 
an unlensed universe with Al — 0.87 ± 1.05. Finally, 
we find A L = 0.96 ± 0.60 and A L = 1.05 ± 0.69 when 
our constraints on the lensed power spectrum using the 
W and V frequency bands respectively are added to the 
WMAP 7-year data and is A L = 0.97 ± 0.47 when the 
W and V frequency bands are combined.. 

The amplitude of the Cf^ reported here is consistent 
with the naive expectation based on a Fisher matrix 
estimate that suggests a measurement no better than 
2.5o\ A variety of effects could be aiding the detection; 
our estimator is sensitive to the full trispectrum while 
a Fisher matrix estimate based on the reduced trispec- 
trum with one pairing could underestimate the expected 
signal-to-noise ratio (see, Fig 3 of Hu 2001b). Effects 
such as the correlation between lensing and secondary 
anisotropies, such as the Sunyaev-Zel'dovich effect, could 
enhance the signal in the lensing trispectrum (Cooray & 
Kesden 2003). While the full WMAP dataset will slightly 



improve the measurement we report here, a confirmation 
of our result showing a detection of the lensing power 
spectrum will come from Planck data, which is expected 
to make a larger than 60cx detection of Al (Hu 2001b). 

As the measurement of Cf becomes more precise, it 
will provide tighter constraints on the cosmological pa- 
rameters than can be obtained using temperature and 
polarization information alone. One obvious future tar- 
get is a measurement of the sum of the neutrino masses, 
leading to a direct cosmological determination of the neu- 
trino mass hierarchy (Kaplinghat et al. 2003; Lesgour- 
gues et al. 2006). A the lensing power spectrum fur- 
ther constrains the dark energy equation of state w as 
well as early dark energy models in general (Kaplinghat 
et al. 2003; Joudaki in Prep.). The ultimate goal of 
CMB polarization measurements is the primordial grav- 
itational wave signal in the so-called B-modes of polar- 
ization (Baumann et al. 2009; Bock et al. 2009). The 
signal however is confused by the lensing effect that con- 
verts a small fraction of dominant scalar polarization in 
E-modes to B-modes (Kesden et al. 2002; Knox & Song 
2002,Zaldarriaga, & Seljak 1998). The fast estimator we 
have presented here for temperature maps can be gen- 
eralized for CMB polarization maps (Munshi et al. in 
prep) and can be used to de-lens CMB B-modes to sep- 
arate the signal from primordial gravitational waves and 
lensing of E-modes. 
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